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We take a closer look at the fundamental Casimir-Polder interaction between quantum particles 
and dispersive dielectric surfaces with surface polariton or plasmon resonances. Linear response 
theory shows that in the near field, van der Waals, regime the free energy shift of a particle contains a 
thermal component that depends exclusively on the population/excitation of the evanescent surface 
polariton/plasmon modes. Our work makes evident the link between particle surface interaction 
and near field thermal emission and demonstrates how this can be used to engineer Casimir-Polder 
forces. We also examine how the exotic effects of surface waves are washed out as the distance 
from the surface increases. In the case of molecules or excited state atoms, far field approximations 
result in a classical dipole-dipole interaction which depends on the surface reflectivity and the 
mean number of photons at the frequency of the atomic/molecular transition. Finally we present 
numerical results for the CP interaction between Cs atoms and various dielectric surfaces with a 
single polariton resonance and discuss the implications of temperature and retardation effects for 
specific spectroscopic experiments. 


The Casimir-Polder (CP) interaction between a polar- 
isable quantum object (atom or molecule) and a surface 
arises from quantum fluctuations in vacuum. It’s an ex¬ 
cellent candidate for fundamental tests of cavity quan¬ 
tum electrodynamics and crucial for any experiments at¬ 
tempting to measure non-Newtonian gravity interactions 
mm ■ CP forces are also relevant in physical chemistry 
playing an important role in the interpretation of physi¬ 
cal phenomena such as atomic adsorption and desorption 
from hot surfaces or even surface chemistry and cataly¬ 
sis . The continuous urge for miniaturisation has led to 
integrated devices, such as atom and molecule chips la¬ 
in., used for a variety of applications and more recently 
tapered nano-fibers were used to trap atoms at distances 
as small as 200 nrn away from the surface CHS, where 
atom-surface forces become exceedingly relevant. Novel 
trapping schemes that exploit the complexity of the van 
der Waals (vdW) potential of excited atoms have also 
been proposed m- 

The most basic description of the CP effect is that 
of a classical dipole interacting with its surface induced 
image. This approach is mostly valid in the vdW ( z _3 
law) regime, but QED theory [TT] revealed that when 
distances are larger than the wavelength corresponding 
to atomic transition, retardation effects scramble the in¬ 
teraction giving a z ~ 4 distance dependence. For ex¬ 
cited state atoms or molecules an additional contribution 
Emus] resembling the interaction of an antenna with its 
own reflected field has to be considered JT4j. Thermal cor¬ 
rections to the CP force are analogous to the black body 
radiation induced corrections to the well known Lamb 
shift In thermal equilibrium the problem has been 
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considered by several authors mm- A novel behaviour 
was predicted when the surface and the vacuum are at 
different temperatures El- 

Accurate experimental demonstrations of the vdW law 
were given in a series of experiments performed with 
beams of Rydberg atoms m as well as spectroscopic 
selective reflection experiments mm ■ Retardation 
effects were also demonstrated (22] with ground state 
sodium beam. Several experiments with cold atomic 
clouds have also been performed mm and a BEC po¬ 
sitioned 6 — 12pan away from a silica surface was used 
to demonstrate the temperature dependence of the atom 
surface interaction out of thermal equilibrium |28j . 

In thermal equilibrium, temperature effects had re¬ 
mained elusive and were only very recently demonstrated 
using spectroscopic measurements in thermal vapour 
cells EE that probe atoms at distances on the order 
of 100 nm away from the surface. Critical to this experi¬ 
ment is the probing of excited state atoms that, much like 
molecules and in contrast to ground state atoms, have 
the advantage of presenting numerous dipole couplings 
in the mid and far infra-red. At these frequencies di¬ 
electrics support surface polariton modes whose thermal 
excitation creates nearly monochromatic electromagnetic 
fields (compared to the well-known black body radiation) 
that evanescently decay away from the surface [ST] . 

Here we use quantum mechanical linear response the¬ 
ory to calculate the thermal CP interaction at all dis¬ 
tances away the surface. By resumming the Matsub- 
ara frequency expansion we derive analytical expressions 
which demonstrate that the thermal component of the 
vdW interaction can be considered as a shift of the atomic 
levels due to near field thermal emission of surface modes. 
As such, its sign and strength depend exclusively on 
the relative position (detuning) of the atomic transitions 
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compared to the frequency of the surface-polariton res¬ 
onance. We show that the thermal excitation of surface 
polariton /plasmon modes can have observable effects 
even for low lying excited state or ground state atoms or 
molecules. We also examine scenarios where temperature 
changes can lead to a complete cancellation or change of 
sign of the atom-surface interaction. This allows engi¬ 
neering of CP forces with the use of temperature. We 
also derive analytical expressions for the resonant contri¬ 
bution to the CP interaction in the far field regime. The 
resonant contribution, dominant in many cases of inter¬ 
est (e.g. excited atoms and molecules) resembles a classi¬ 
cal dipole-dipole interaction and depends on the surface 
reflectivity. Unlike the non-resonant CP predicted be¬ 
haviour for ground state atoms [32] , the resonant contri¬ 
bution preserves a strong anisotropic component, charac¬ 
teristic of dipole-dipole interactions. We finally discuss 
the transition between the two regimes using numerical 
calculations. We demonstrate that controlling the CP 
interaction with temperature is no longer possible as the 
distance from the surface increases. We also show that 
retardation effects are relevant even for spectroscopic ex¬ 
periments at nanometric distances away from the surface. 


I. CASIMIR-POLDER INTERACTION 

We start by considering the CP free energy shift of a 
quantum particle A F a at a given energy level |a), which 
is the sum of a resonant and a non-resonant contribution 
EH 133]. We follow the formalism of M-P Gorza et. al. 
m, describing the free energy shift at a finite tempera¬ 
ture T. 


AF 0 = AF; + AF a " 


(1) 


To simplify the notation we will expand our reasoning 
for a two level |a), \ b) system. For a real multilevel sys¬ 
tem one simply has to sum all the contributions of all 
individual dipole couplings. 

At a finite temperature the non resonant term is given 
by the following sum : 


Ar r = 


2 ksT 


OO f 


E (*aHFG afi (z,iS k ) 

k =0 



( 2 ) 


case, the linear susceptibility gives the reflected displace¬ 
ment field at a point r due to a dipole /2(w), oscillat¬ 
ing at a frequency w, positioned at r', via the relation 
—* —> ^ ^ 

D(r,r',u)) =G (r,r',uj)p,(uj). In our case G is evaluated 

for f = r', because we’re interested in dipole interacting 
with their own reflected field. Due to the cylindrical sym- 

metry G is only a function of frequency and the distance 
z of the dipole from the reflecting wall. 

The resonant part of the CP shift is written as HU: 

A K = n(u 0 ,T)^nfRe [G aP (z, \w a \)] (3) 

The mean occupation number n(u) 0 ,T) = —- of a 

e k B T -1 

mode according to Bose-Einstein statistics is here ex¬ 
tended to negative frequencies. Note that in the case 
of virtual emission uj 0 < 0 the sign changes and an ad¬ 
ditional photon due to spontaneous emission is added, 
n(w 0 , T) = - [1 + n(-w 0 , T)}. 

The problem reduces to calculating the linear suscep¬ 
tibility function. We will focus on the simple geometry 
of a semi-infinite surface with a bulk dielectric constant 
e(w) and an atom that is in the vacuum . The diagonal 
terms of the linear susceptibility [HUM calculated 
for an imaginary frequency are : 

c3 rco 2 , iu 

G xx (z, i£) = G yy (z, *0 = ^ j « ( u 2 R p - R s ) du 

(4) 

£3 /»oo 

G zz {z,i£) = J (u 2 -l) RPdu (5) 

In the above equations R p and R s are the Fresnel re¬ 
flection coefficients and u is a dummy integration vari¬ 
able. The problem has no simple analytical solutions 
apart from the famous case of an ideal conductor or an 
ideal dispersion-less dielectric. For a real surface one has 
to resort to numerical simulations except in the limiting 
cases when 1, i.e the long range case where re¬ 

tardation effects are important, or 2 ^ z <C 1 which is the 
electrostatic or van der Waals regime. 


II. VAN DER WAALS INTERACTION 


We use the Einstein notation , implying a summation 
over the index variables a and /? that denote the Carte¬ 
sian coordinate components. The prime symbol signi¬ 
fies that the first term of the sum should be multiplied 
by 1/2. The transition frequency oj 0 = (Et, — E a )/h 
depends on the energy difference between the two lev¬ 
els. It takes positive or negative signs depending on the 
nature of the coupling (absorption or emission). Also 
= 27 t ^|—k are the Matsubara frequencies, and 
are the dipole moment matrix elements and G a p(z,it;k) 
are the components of the linear susceptibility matrix 
of the reflected field, defined in mis]- In the general 


Since the first experimental demonstration of retarda¬ 
tion effects [20] the vdW interaction has been mostly con¬ 
sidered as an electrostatic limit of the CP interaction. A 
QED description is however necessary when the atom is 
in the presence of a hot surface. Here we pursue further 
the results of m by considering real dielectrics with one 
polariton resonace in order to illuminate the physics of 
CP interactions in the presence of polaritons or plasmons. 
We show that excited surface waves create intense ther¬ 
mal fields in the vicinity of the surface (near field thermal 
emission) that offer a way to control the vdW interaction, 
in a way that is not possible in the far-field regime. 
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In the near field, vdW regime (z < jj, where A 0 is 
the atomic transition wavelength) the linear susceptibil¬ 
ity has an analytical solution given by pT5i : 


— Gyy(z,Uj) — 


G zz (z,u) 


2 

( 2 # 


(2z) 3 5(w) 

(6) 

S{u) 

(7) 


Eqns. (6,7) are valid for both real and imaginary frequen¬ 
cies when 2 ' u, ' z <C 1. Here S is the frequency dependent 
image coefficient which is a function of the bulk dielectric 
constant e(ui). It is given by: 


5V) = 


eM ^ 1 

e(w) + 1 


( 8 ) 


We start by using a single resonance model to describe 
the dielectric constant of the surface. 


e (w) = e inf + 


(^st 

— UJ 2 — iTui 


(9) 


where ei n f and e s t are constants giving the dielectric 
constant at the two extreme ends of the spectrum, T 
is phenomenological decay rate and wj- is the transverse 
frequency of oscillations. The above equation models di¬ 
electrics with one surface polariton resonance. The valid¬ 
ity of eqn. (9) is limited to a certain frequency range (see 
also discussion in [IE]), however, it accurately describes 
the CP interaction between many atom-surface systems. 
A more realistic model should account for multiple reso¬ 
nances, however this scenario will not be considered here. 

We first assume that T = 0. In this case, the surface 

polariton frequency is given by los = wt e^f+i anc ^ 
image coefficient becomes: 

5 ( W ) = s inf + (s st - s inf ) “ s 2 (io) 

UJ g — UJ 


Here S ln f and S s t represent the values of S(ui) calculated 
for a n f and e st respectively. Eqn.pT| diverges at the 
polariton frequency but this is a small price to pay for 
keeping our analytical expressions simple with a clear 
physical interpretation. 

Reporting the linear susceptibility eqns. j 6 j [7| and 
the dielectric constant eqn. § into eqn. ([ 2 ]), we can 
sum the Matsubara frequency expansion and arrive to 
the following analytical expression for the non-resonant 
free energy shift of the vdW interaction: 


the relation coth ^ 2 ^ T ^ = 2 n(co,T) + 1. One notices 

the characteristic z~ 3 dependence of the vdW interac¬ 
tion. The terms inside the brackets in eqn. (Ill depend 


only on the dielectric properties of the surface and on 
temperature. We refer to them as the non-resonant part 
of the image coefficient. The seemingly complicated tem¬ 
perature dependence of the non-resonant term, given by 
a complex sum over imaginary frequencies, essentially re¬ 
duces down to the number of photons in the atomic and 
polariton frequencies. On the other hand the resonant 
part of the shift, eqn. ([3|, only depends on the number 
of photons in the atomic frequency and is given by: 


a f: = - 


Gl c 


Re [5(w 0 )] — Re [jS(w d )] coth 


/ hu 0 \ 


\2k B T J 

( 12 ) 

We refer to the term inside the brackets in eqn. © as 
the resonant part of the image coefficient. A simple in¬ 
spection of the above equations reveals that the temper¬ 
ature dependence of the resonant contributions cancels 
out with the first term of the non-resonant contribution. 
Adding eqn. © and eqn. © we find the total free 

energy shift A F a = -%-r (co a , T) = Here, C 3 is 

the vdW coefficient and r (w OJ T) is the image coefficient 
given by: 


r (uj 0 ,T) = Re [5(w 0 )]-l- 


(Rst ^inf) 


coth 


W- - w 


/ hu s \ 
\2k B T ) 
(13) 

The above equation gives the temperature dependence 
of the vdW interaction. It depends strictly on thermal 
fluctuations at the polariton frequency i.e on the ther¬ 
mal excitation of the evanescent surface polariton modes 
and not on the number of thermal photons at the atomic 
frequency. In the two extremes of the spectrum, when 
uj 0 los and lo 0 <C los all temperature dependence 
vanishes, in accordance with previous results given for 
an perfect conductor j.T 8 j. In consistence with the clas¬ 
sical picture, when material dispersion is neglected the 
vdW attraction is independent of temperature. Similar 
results can be obtained with the methodology of m but 
eqn. (13) is more general since it is valid for all temper¬ 
atures, including T=0 and includes both virtual absorp¬ 
tion (w 0 > 0) and emission (uj 0 < 0). In the limit of high 
photon number the temperature dependent part of eqn. 
(13) can also be derived by calculating the Stark shift 
induced by the thermally populated evanescent surface 
polariton modes whose density of states, p B , is given by 

EH E2: 


P C 

A F2 r = - URe [SM coth 


/ hup \ 

\2 k B Tj 


( Sst Sj n j ) LOpUJs 


coth (AM] 

\2.k B T J 


( 11 ) 


Here the constant Cf c = represents the vdW 

coefficient for a perfect conductor. We have also used 


Pd{v,T,z) 


1 

87 t 2 loz 3 


Im 


e-l 
£ + 1 


(14) 


This last approach is analogous to the one used in [3S] 
in order to calculate the atomic shifts due to black body 
radiation. 

The calculation becomes significantly more cumber¬ 
some if T 7 ^ 0 , but analytical expressions can also be 
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FIG. 1. (Color online)(a) The real part of the surface re¬ 
sponse for real frequencies. The results of the hydrodynamic 
model,eqn. Q are shown for sapphire (red or dark grey) and 
BaF 2 (grey). The black dotted (sapphire) and dashed (BaF 2 ) 
lines represent the results of the elaborate model of ref. [36]. 
(b) The surface response for imaginary frequencies using the 
exact same colour/line coding. The dielectric constant and 
therefore the surface response are strictly real in this case. 


found. In most cases of practical interest (r -C ws) 
many terms can be neglected and the image coefficient 
becomes: 


r (w OJ T) = Re [S(w 0 )] + 

2W„WS {Sat - Sinf ) [r 2 + 2(w 2 - w|)] / hu S 

[r 2 + 2(w 2 - w|)] 2 + 4r 2 w| co V 2 k B T 

(15) 


Within the limits of eqn. (15) the cancellations that make 
the vdW interaction independent of the number of pho¬ 
tons at the transition frequency are still valid. The 
above approximation deviates from the exact solution in 
the vicinity of the polariton resonance only by a few per¬ 
cent, whereas in the rest of the spectrum differences are 
negligible. 

To illustrate the importance of these results we start 


by giving in Table I a list of all the parameters involved 
in eqn. (|9| i.e e st , ei n f , cot and T for some dielectrics. 
Amongst them, sapphire is very commonly used in atom- 
surface interaction experiments whereas CaF 2 and BaF 2 
have been considered for this purpose (see [ 55 ]| 55 ]) due 
to their isolated surface resonances in the relatively far 
infra-red. SiC is most commonly used for near field ther¬ 
mal emission measurements mm- Here we explore the 
potential interest of performing CP measurements with 
this material. We also restrict ourselves to dielectrics 
whose dielectric constant is adequately described by the 
hydrodynamic model throughout the visible, near and 
far infra-red frequency range. In the case of CaF 2 , BaF 2 
and sapphire the parameters are deduced by fitting eqns. 
( |8|9[ ) to the experimental data given in ref. [ 35 ] , whereas 
for SiC the parameters are taken directly from [31 j . In 
Fig0 a) we plot the real part of the surface response 
vs real frequencies (w/27t) for both sapphire and BaF 2 . 
The results of eqn. Q (coloured lines) are compared 
to those of the a more elaborate model described in [ 35 ! 
(black lines). The surface response for imaginary fre¬ 
quencies (£/27r) is shown in Figjljb). The hydrodynamic 
model of eqn. © reproduces very well the surface re¬ 
sponse of the above dielectrics. The only observable 
discrepancy between the two models is due to the ex¬ 
istence of an additional, albeit much smaller, surface res¬ 
onance at 480 cm -1 , in the case of sapphire. One should 
not be deceived into thinking that eqn. ([9]) is a perfect 
model of the dielectric constant itself. However, Figjl] 
clearly demonstrates that our simple model can be reli¬ 
ably used to predict the van der Waals free energy shifts 
of atoms/molecules against any of the dielectric surfaces 
of Table I. 


TABLE I. Parameters used to model the dielectric constant 
of sapphire, CaF 2 , BaF 2 and SiC. 



€inf tst 

If (cm -1 ) 

If (cm -1 ) 

2tiT 

UJS 

BaF2 

2.12 7.16 

291 

179.9 

0.047 

CaF 2 

2.02 6.82 

416.2 

258.6 

0.063 

Sapphire 3.03 9.32 

828.9 

518 

0.02 

SiC 

6.7 10 

947.8 

793 

0.005 


Having established the validity and the main features 
of our model we can proceed to calculate the interac¬ 
tion between some realistic atom surface systems. In 
'Fig. m we plot the C 3 coefficient as function of tem¬ 
perature for Cs*(7D 3 / 2 ) against sapphire and SiC. The 
C 3 coefficient for this atomic level depends strongly on 
the 7D 3 / 2 —► 5 F 5/2 dipole coupling at 923cm -1 . The 
sapphire resonance at 829cm -1 is on the red side of this 
transition frequency leading to an increase of the vdW co¬ 
efficient as a function of temperature as demonstrated ex¬ 
perimentally in ref. [30]. Conversely, the same atom near 
a SiC surface displays a completely different behaviour. 
The SiC resonance at 945cm -1 , is on the blue side of 
the transition frequency, thus the vdW coefficient now 
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FIG. 2. (Color online) The vdW coefficient as a function of 
temperature for a Cs(7-D 3 / 2 ) atom against a sapphire (red or 
dark grey) and a SiC (black) surface. 


decreases with temperature. According to our theoret¬ 
ical estimates the fundamental long-range atom-surface 
interaction should be null at T ~ 1200A', a temperature 
range which could be experimentally achievable. Fig. [2] 
shows the possibility of controlling the CP interactions 
close to dispersive surfaces via temperature. Other atom- 
surface systems have been considered for achieving repul¬ 
sive vdW potentials j5B] at finite temperatures but so far 
there has been no experimental proof of this effect [29] . 

It is also of interest to estimate the effects of surface 
polaritons for ground state atoms. This situation is rel¬ 
evant for most experimental measurements using with 
cold atoms [2] [23J [20] US] ■ The finite temperature correc¬ 
tions (T = 300AT) in the near field regime are in this case 
0.02%, 0.16%,1.29% and 1.75% for SiC, sapphire ,CaF 2 
and BaF 2 respectively. These numbers indicate that near 
field thermal emission has negligible effects compared to 
the experimental precision of most experiments so far, 
but it could have implications in the case of precision ex¬ 
periments aiming at putting new limits to the existence 
of Non-Newtonian gravity forces mm- The above num¬ 
bers should be considered as indicative because we have 
used the simplified model of eqn. ®> and they are calcu¬ 
lated in the vdW electrostatic limit , which restricts their 
validity to very small distances away from the surface. 


III. FAR FIELD APPROXIMATION 

The behaviour of the non-resonant term of eqn.pj) has 
been discussed extensively in the past (see e.g. [18ll28| f39 ] 
and references therein). We’ll just remind here that at 
T = 0 there’s a passage from a z~ 3 to a z~ 4 law for 
the energy shift at large distances from the surface. At 
a finite temperature T there’s a second cross over to a 
z~ 3 law when z , where A t = is usually 

referred to as the thermal wavelength. In this case the 


sum of eqn. ([2]) is dominated by the first term |39| and 
therefore gives: 


A FT 


2 e(0)-l k B T Cl c 

e(0) + 1 huj 0 z 3 


(16) 


In the far field the free energy shift is independent of 
material dispersion. Only on the DC or ’static’ value 
of the dielectric constant at zero frequency come into 
play. This far field or high temperature limit is sometimes 
refered to as the Lifshitz regime [28]. 

The behaviour of the resonant term in the far field is 
rather different. Using the imaginary frequency formulas 
above we can see that in the far field the exponential de¬ 
cays rapidly. The multiplying functions ( u 2 R p — A s ) and 
( u 2 — l) R p vary rather slowly and one can find approx¬ 
imate expressions for the linear susceptibility by Taylor 
expanding them around 1 to the lowest non-zero order. 
The real part of the linear susceptibility for real frequen¬ 
cies w 0 is the relevant quantity for the resonant term of 
the free energy shift, which is given by: 


Ae \Gxx,yy (-A )] — n |p(^o) \cosi^2h 0 z 0 -f- t^>(cu 0 )) (17) 

Zz n 


Re[G zz (z, w 0 )] = -^\p(ui 0 )\sin(2k 0 z 0 + 4>(u> 0 )) (18) 

where p(w 0 ) = \p{(jj 0 )\e % ^ u °> is the frequency dependent 
complex reflectio n coeffi cient of the surface at normal 
incidence. Eqns. (17 18) are valid when z where 

A 0 is the transition wavelength. As has been pointed out 
before M the resonant term resembles the interaction of 
an antenna with its own reflected field, which oscillates 
between repulsion and attraction with a period of A 0 /2. 


In the case of a perfect conductor, eqns.(17 18) lead to 
the cavity QED shifts reported in P3]. The term is highly 
anisotropic since a dipole antenna does not radiate on its 
axis. For this reason Re[G zz (z,uj 0 )\ decays much faster 
than Re[G XXtVy (z , w D )]. In reality this term will dominate 
the far field of CP interaction for excited state atoms 
or molecules 00 SB. At very large temperatures or at 
very large distances from the surface it will eventually 
dominate the interaction between a surface and a ground 
state atom. 

In Fig[3]we plot the reflection coefficient as a function 
of frequency. The limitations of the hydrodynamic model 
are here visible, especially for sapphire. The discrepan¬ 
cies are due to additional resonances that have been ne¬ 
glected 021. The most noteworthy feature of Figgis that 
the reflection coefficient does not display the same varia¬ 
tions as the surface response which changes sign around 
the surface polariton frequency (Fig. [l]). This shows that 
exotic effects such as resonant vdW repulsion 03] or tem¬ 
perature effects due to the thermal excitation of surface 
polaritons [30] should be washed out as the distance form 
the surface increases. 
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FIG. 3. (Color online)The absolute value of reflection coef¬ 
ficient as a function of frequency. The solid lines represent 
the reflection coefficient as calculated using the hydronamic 
model of eqn. |9j in the case of sapphire (red or dark grey) 
and BaF 2 (grey). The black dotted (sapphire) and dashed 
(BaF 2 ) lines represent the results of the model described in 
ref. j|g. 


IV. DISCUSSION 


In order to demonstrate the effects of retardation and 
thermal excitation of the surface we focus our attention 
in some concrete examples. We first choose two cou¬ 
plings (virtual emission) with frequencies —820c?n -1 and 
—838cm -1 corresponding at the two extrema of the sap¬ 
phire image coefficient, close to the sapphire resonance at 
829cm -1 (Fig(T|). In Fig. [4^, we plot the normalised CP 
shift due to these couplings as a function of distance from 
a sapphire surface. The shift is calculated by numeri¬ 
cally integrating eqns. @§ using the hydrodynamic 
model for sapphire’s dielectric constant, without any ap¬ 
proximations. In the near field the two dipole couplings 
present shifts of opposite sign. The first (blue line at 
—820cm -1 ) corresponds to a large attraction whereas the 
second (red line at —838cm -1 ) corresponds to a repul¬ 
sion. This exotic effect persists only in the nanometric 
scale, at distances smaller than 1/zm. The two curves 
converge for larger distances, oscillating from attraction 
to repulsion in a very similar fashion. The reader should 
note that even though the amplitude of these oscillations 
is almost the same, there is a noticeable phase difference 
between them. The estimated free energy shifts using 
the far field approximation of eqns. (17 18) are also 


shown as black dashed and dotted lines. At large sep¬ 
arations (z > 10/rm)tliey coincide almost perfectly with 
the numerical calculations. In Fig. & we plot the free 
energy shift due downward coupling at — 830cm -1 for a 
BaF 2 surface at two different temperatures T = 200 K 
and T = 600A7 Due to the thermal excitation of the 
BaF 2 surface polariton at 291cm -1 the vdW interaction 
changes sign from attraction to repulsion. The near field 



distance (pm) 


FIG. 4. (Color online)(a) The free energy shift multi¬ 
plied by z 2 for —820cm -1 (blue or black) and — 838cm -1 
(red or dark grey) downward couplings against a sapphire 
surface at T=300 K. The dashed and dotted lines represent 
the fre e energ y shift as given by the far field approximation 
of eqn (17 181 (b) The free energy shift multiplied by z 2 for 
a — 830cm 1 downward coupling against a BaF 2 surface for 
T = 200A' (blue or black) and T = 600A' (red or dark grey) 
lines. 


temperature dependence is governed by the number of 
thermal photons at the the polariton frequency ujs of 
BaF 2 , as described by eqn. (13). Conversely the temper¬ 
ature dependence in the far field is related to the number 
of thermal photons at the transition frequency u 0 , as can 
be seen by eqns. ( |3|17|18D . Here again the effects of the 
evanescent polariton modes are only present in the nano¬ 
metric range|44? ]. The transition frequencies used in 
Fig|4] are in the vicinity of the 7P —► 6D transitions of 
both Cs and Rb. The vdW interaction of Cs*(6D 3 / 2 ), 
Rb*(6D 3 / 2 ) and Rb*(6D 5 / 2 ) was experimentally investi¬ 
gated in the past with selective reflection experiments in 
vapour cells j45] . 


Finally, we perform a complete calculation of the 
Casimir-Polder free energy shift for the case of ground 
state Cs(6S!/ 2 ) and the low lying excited state 
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TABLE II. Individual contributions to the C 3 coefficient 
of the most important dipole couplings for Cs( 6 Si/ 2 ) and 
Cs( 6 Pi/ 2 ) close to a perfect conductor. Transition probabili¬ 
ties are taken from ref. 46j 

Cs( 6 S!/ 2 ) A(/tm) w(cm _1 ) C^ c (kHz/j,m 3 ) 

6 P i/2 0.894 11178.24 0.94 

6 P 3/2 0-852 11732.35 1.62 


Cs(6P 1/2 

) A (/wre) 

cj(cm J ) 

C^ c (kHzfim,' 

6 S 1/2 

-0.894 

-11178.24 

0.94 

5D 3 /2 

3.01 

3321.25 

1.64 

7S 1/2 

1.36 

7357.27 

0.59 

6 D 3/2 

0.876 

11410.65 

0.64 

7D 3 /2 

0.673 

14869,62 

0.15 


Cs*( 6P 1/2 ). Here we take into account only the most 
important dipole couplings (Table II), ignoring the con¬ 
tribution of core excitations |4T]. In Fig { 5 ] we show our 
results, focusing first (Fig{5^t) on rather short distances 
from the sapphire wall. There is a dependence of the 
CP interaction on temperature, which is more impor¬ 
tant for the Cs*( 6P!/ 2 ) state because the 6P!/ 2 -t 5-D 3 / 2 
coupling is closer to the sapphire polariton frequency. 
The effect should be barely measurable with the preci¬ 
sion of spectroscopic atom-surface interaction [48] exper¬ 
iments. One notices that for the ground state Cs atom, 
the van der Waals limit is not yet reached even at such 
short distances, characteristic of the distance dependence 
of the non-resonant contribution to the CP interaction 
mmmm- Conclusive experimental evidence of this 
behaviour are only given for distances greater than 150 
nm [25 . This is less true for the Cs*(6Pi/ 2 ) which seems 
to follow a z~ 3 dependence for distances smaller than 
lOOnm. Intuitively one can assume that this is because 
the wavelength of the dominant coupling is larger for the 
6 Pi/ 2 state (3/tm instead of 0.852/im for 6 Si/ 2 ), but this 
is not true. Remarkably it is the interplay between the 
different couplings (some virtual emissions and some vir¬ 
tual absorptions) that leads to this phenomenological ad¬ 
herence to the vdW z~ 3 law. This coincidence could have 
important implications for spectroscopic experiments in 
the near held of the atom surface interaction, which are 
sensitive to the free energy difference between levels. For 
low lying excited states (such as the 6 Pi/ 2 ) experimen¬ 
tal measurements of the CP shift could vary depending 
on the typical probing distance. In this particular ex¬ 
ample the 6 Si/ 2 —> 6Pi/ 2 shift is almost doubled be¬ 
tween 0 and 100 nm and is clearly not following the vdW 
z~ 3 law. This should be taken into consideration for ex¬ 
periments on low lying excited states of alkalis pTST(5(7l 
(see also Discussion and Prospects section of [51]). It’s 
worth mentioning that in the selective reflection experi¬ 
ments reported in [48] the measured C 3 coefficient is 1.4 
kHz /mi 3 , a value which is more consistent with Fig. [5] 




FIG. 5. (Color online)(a) The CP free energy shift multiplied 
by z 3 for Cs( 6 S 1 / 2)5 Cs*( 6 Pi/ 2 ) against a sapphire surface and 
the difference between them at two different temperatures, 
T = 500R' (black) and T = 1200 K (red or dark grey) for 
distances between 0 and 200 nm. (b) The CP free energy 
shift multiplied by z 3 for Cs( 6 Si/ 2 ) against a sapphire surface 
for a distance range between 0 and 20 /tm, at two different 
temperatures, T = 500 K (black) and T = 1200A' (red or dark 
grey). The dashed and dotted lines show the results of the 
far field approximation of eqn. (16). For higher temperatures 
eqn. (16) is valid for distances closer to the surface. 


(for distances greater than 50-100 nm) than the theo¬ 
retical vdW prediction of 0.8 kHz /mi 3 (essentially for 
z=0). Conversely high lying excited states exhibit huge 
CP shifts compared to low lying states and should be 
less affected by this problem [22] 2D, 30] FT5l ? ]. In Fig. 
Ep we plot the CP shift of the Cs ground state atom at 
a much larger range of distances (0-20/im). The non¬ 
resonant term dominates and the approximation of eqn. 
|16| is mostly valid in the far field. At very high temper¬ 
atures (T=1200 K) one begins to see the beat between 
the QED oscillations of the resonant term at 894nm and 
852nm. Putting these predictions to a real experimental 
test is at the moment extremely challenging. 




























V. CONCLUSIONS 


We have analysed the thermal effects of the CP inter¬ 
action when the surface is at thermal equilibrium with 
the surrounding environment. We derived simple analyt¬ 
ical expressions in the case of a dielectric with one surface 
resonance that provide a transparent physical interpreta¬ 
tion of the temperature dependence of the CP free energy 
shift of atoms or molecules. Our work shows that in the 
near field thermal effects are entirely due to the exci¬ 
tation of evanescent surface modes (near field thermal 
emission). The results can be easily extended to the case 
of metals where plasmon frequencies are typically at UV 
wavelengths. Using more realistic models for dielectrics 
is also straightforward but significantly more tedious. We 
also derive simple expressions in the far field approxima¬ 
tion, valid not only for ground state but also excited state 
atoms and molecules. We show that retardation effects 
can be significant for spectroscopic experiments on low 
lying states performed at nanometric distances away from 
the surface. This raises important questions concerning 
the validity of the z~ 3 vdW law. It is often neglected that 
the CP interaction can also affect the radiative proper¬ 
ties (transition rates) of atoms or molecules [MIES]. Our 


conclusions can be extended to calculate a distance de¬ 
pendent transition linewidth da hi using the imaginary 
part of the linear susceptibility matrix. Although our 
treatment was entirely performed in the case of thermal 
equilibrium, our demonstration that the near field tem¬ 
perature dependence is solely due to the excitation of 
evanescent surface waves suggests that in the near field 
our results are also valid in an out of equilibrium case 
with T being the temperature of the surface [H] . A com¬ 
plete treatment of the out of equilibrium case for excited 
state atoms or molecules is, however, more challenging. 
Finally, the simplicity of our results could render them 
very useful when applied to intrinsically more compli¬ 
cated problems such as multi-layered dielectrics [13], or 
2-D dielectrics such as graphene deposited on dielectric 
substrates [53] or even meta-surfaces. 
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